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Abstract —Wireless underground sensor networks (WUSNs) 
present a variety of new research challenges. Magnetic induction 
(MI) based transmission has been proposed to overcome the very 
harsh propagation conditions in underground communications 
in recent years. In this approach, induction coils are utilized 
as antennas in the sensor nodes. This solution achieves longer 
transmission ranges compared to the traditional electromagnetic 
(EM) waves based approach. Furthermore, a passive relaying 
technique has been proposed in the literature where additional 
resonant circuits are deployed between the nodes. However, 
this solution is shown to provide only a limited performance 
improvement under practical system design contraints. In this 
work, the potential of an active relay device is investigated which 
may improve the performance of the system by combining the 
benefits of the traditional wireless relaying and the MI based 
signal transmission. 


I. Introduction 

The objective of Wireless Underground Sensor Networks 
(WUSNs) is to establish an efficient wireless communication 
in the underground medium. Typical applications for such 
networks include soil condition monitoring, earthquake predic¬ 
tion, border patrol, etc. id, ia. Since the propagation medium 
is soil, rock, and sand, traditional wireless signal propagation 
techniques using electromagnetic (EM) waves can be only 
applied for very small transmission ranges due to a high path 
loss and vulnerability to changes of soil properties, such as 
moisture or water contents 0. 

Magnetic induction (Ml)-based WUSNs were first introduced 
in El, 0] and make use of magnetic antennas implemented 
as coils. In recent years, some system models have been 
developed in order to investigate the system behavior and 
typical properties of the MI based transmission channels, e.g. 
0 , 0 , and 0 . More realistic channel and noise models for a 
point-to-point transmission were derived in ||7j. These models 
incorporate the losses due to the transmission medium and 
the power reflections between the coils. Furthermore, it was 
shown that the transmission through a conductive medium like 
soil can only be established using a carefully optimized set 
of system parameters. Recently, novel modulation approaches 
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for MI based transmissions using practical signal processing 
components (transmit, receive, and equalization filters) have 
been proposed in (8). 

The idea of deploying additional resonance circuits (passive 
relay devices) between two transceiver nodes and combining 
them into waveguide structures has been first proposed in 
m for MI based communications systems and in m for 
MI-WUSNs. Similar to traditional wireless relaying concepts 
this approach is supposed to benefit from a lower path loss. 
However, as it was pointed out in 0, the MI waveguides 
show limited performance in the underground medium under 
realistic design constraints due to a very narrow bandwidth, 
as discussed in detail in J8]. Hence, in this work, we propose 
to add a further functionality to the relay devices, such that 
the received signal at the relay is actively processed and 
retransmitted. This additional functionality for MI enabled 
communication systems is already mentioned in DSD- How¬ 
ever, the potential of this technique for the MI-WUSNs has not 
been shown and analyzed as of to date. Hence, we investigate 
the properties of active relaying schemes in MI-WUSNs. 
Furthermore, we consider the possibility of establishing a full 
duplex relaying mode, which may dramatically improve the 
performance of this scheme. 

This paper is organized as follows. In Section [II] the system 
model is derived, which incorporates the influence of a single 
active relay onto the signal transmission and reception. To 
this end, we investigate the frequency-selective path loss of 
the relayed transmission and the colored noise at the receiver 
front-end. In Section III the problem of system parameter 
selection and power allocation is considered. Numerical results 
are provided in Section [IV] and Section [V] concludes the paper. 


II. System Model 

In this work, we derive new channel and noise models for 
direct MI transmission in presence of a single active relay. We 
assume a simple MI relay network with two MI transceivers 
(called source and destination) and one relay between them. 
All three devices contain the same set of passive elements of 
the resonance circuit: 



receive power spectral density is given by 



Fig. 1. Equivalent circuit model for the source-relay-destination link. 


• an induction coil with inductivity L, which is modeled 
as a multilayer air core coil, 

• a resistor with resistance R, which models the parasitic 
wire resistance, 

• a capacitor with capacitance C, which is tuned to make 
the circuit resonant at the carrier frequency fo = 

• a load resistor Rl, which is optimized for minimum 
power reflections at the transceivers for transmissions at 
resonance frequency. 

The basic structure of the relay network is shown in Fig. [T| 
The coupling between the source and the relay is determined 
by the mutual inductance M$r, and between the relay and 
the destination by the mutual inductance Mrd, respectively. 
These mutual inductances are calculated similarly to 0, 

Msr = N 2 3 JsrGsr, ( 1 ) 

ir SR 
a 4 

Mr D = fiTrN 2 — 2 — JrdGrd, (2) 

d n 


where rji denotes the distance between the specified two coils, 
a stands for the coil radius, N is the number of windings, and 
// denotes the permeability of the medium. represents the 
polarization factor according to CD. Gj ) is an additional loss 
factor due to eddy currents as pointed out in 0. We assume 
that all devices are deployed in a homogeneous conductive 
environment (soil) with constant properties over space and 
time. 

Similarly to the previous work (cf. 0), we obtain the transmit 
power spectral density at the source node 

Z1 + (27T fM RD ) 2 
Z g (Z^ + (2nf)^M 2 R + M 2 RD )) ’ 


Pst(f) = 




where Z g = j2nfL + j 2 ^fc + R + Rl denotes the sum of all 
circuit impedances and Ug is the input voltage at the source. 
Correspondingly, Ur stands for the input voltage at the relay. 
The receive power spectral density at the relay is given by 
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Hence, the channel power gain can be calculated as 
\HsR(f)\ 2 = Pgt(f) us i n g 0 anc * 0- addition, the 
transmitted signal from the source is received at the destination 
via passive relaying, which has been investigated in the context 
of MI waveguides in the previous works. The corresponding 
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such that the channel power gain from source to destination 
can be given by \Hgn } p(/) | 2 = When the relay 

transmits its signal, we obtain for the transmit power spectral 
density at the relay and for the receive power spectral density 
at the destination, respectively. 
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This yields a relay to destination channel power gain 
|iW/)| 2 = ^ A major difference compared to the 
relaying in traditional communication systems is an additional 
path loss, which basically results from the linear mapping of 
the received signal at the load impedance onto the transmit 
signal Ur. This mapping can be described using an additional 
channel gain 


\H RR (f )\ 2 


Rl, 
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such that for nonregenerative relaying schemes (e.g. amplify- 
and-forward relaying) the total channel power gain of the 
active relay link (without the amplification A of the relay) 
can be given by 


\H S D, a (f)\ 2 = \H S R(f)\ 2 \H RR (f)\ 2 \H RD (f)\ 2 - (9) 


For the noise modeling, we focus on the thermal noise pro¬ 
duced by the resistors in all three resonance circuits (voltage 
sources Un,s> Un,r. and Un,d in Fig- [TJ- The noise power 
spectral density at the relay and destination, respectively, can 
be calculated similarly to 0 
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where K « 1.38 • 10 -23 J/K is the Boltzmann constant and T 
is the temperature in Kelvin (T = 290 K in this work). 

Since our main focus is on the information theoretic analysis 
of the proposed scheme, the performance metrics are related to 
the achievable data rates, which are calculated using Shannon’s 
channel capacity formula. 











































III. Active Relaying 


o 


In this section we investigate the properties of different 
active relaying schemes such as amplify-and-forward (AF), 
frequency-selective amplify-and-forward (also called filter- 
and-forward, FF) or decode-and-forward (DF) known from 
the literature (e.g. Ifl2l . ifTSl ). In the commonly used half 
duplex relaying mode m, the transmission and reception data 
streams are separated e.g. by means of time-division duplex 
(TDD). However, the corresponding loss in data rate can be 
reduced by applying a full duplex approach. For simplicity, 
we restrict ourselves to the half duplex (HD) relaying mode 
in Section |HI-A| - Section |III-C| In Section |HI-D| the modifi¬ 
cations of the system are discussed, which are needed in order 
to enable the full duplex (FD) relaying mode. 

There are several system parameters, which can be optimized 
for any relaying scheme. Some parameters are typical for the 
MI based transmissions, as discussed in 13, e.g., the resonance 
frequency f 0 and the number of coil windings N. Other 
parameters are typical for the relayed transmission, e.g. the 
transmit power at the relay Pr or the relay position relatively 
to the transceiver nodes. All these parameters (except for the 
design of the transmit filters) are optimized via a full search 
similarly to J3- The optimization of the transmit filters for 
different relaying schemes is discussed below. 


A. Amplify-and-fom’ard relaying (AF) 

The simplest relaying scheme is AF relaying. Here, the 
signal received in the first time slot at the relay is amplified 
using a frequency flat constant A and sent to the destination 
in the second time slot. A is related to the available relay 
transmit power Pr and can be given by 
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where B represents the chosen signal bandwidth. The total 
received power spectral density at the destination in the second 
time slot is 


PdMI) = A 2 ( Pstif ) \H S R(f)\ 2 + PRn(fj) (13) 

X \H RR (f)\ 2 \H RD (f)\ 2 + P D n(f) 

= A 2 P S t(f)\H SD , a (f )\ 2 
+ A 2 \H RR (f)\ 2 \H RD (f)\ 2 P Rn (f)+P Dn (f). 


In order to improve the performance of this scheme, we 
additionally employ the signal received from the source at the 
destination in the first time slot via passive relaying, which 
has a power spectral density 

P Dr i(f) = Pst(f) \Hsd, p (/)| 2 + PduU)- (14) 

The received signals from two time slots may be affected 
by different channel transfer function. Therefore, we combine 
them coherently using the optimal approach of maximum ratio 
combining (MRC). For this, we employ the respective matched 
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Fig. 2. Normalized SNR MRc(f) of the received baseband signal at the 
Destination using AF relaying for different positions of the Relay. 


filters at the destination. The data rate of the relay network can 
be given by 


Caf,hd = 7 ; J 1 °S 2 (1 + SNRmac(/)) d/, (15) 

where we define SNR MRc(f) = SNRi(/) + SNR 2 (/) with 
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Since Pst(f) can only be optimized for a fixed value of A 
and A depends on Pst(f) according to ( | 1 2j i, we propose an 
iterative algorithm in order to maximize the achievable data 
rate. Starting with a uniform power distribution, the value of A 
is calculated in each iteration based on ( [T2| ). Then, the optimal 
power allocation Pst(f) is calculated using the waterfilling ap¬ 
proach. Since SNRmac(/) has a shape Pst(f)K(f), where 
I\(f) is independent of Pst(f ), the waterfilling solution can 
be given by 

Pstif) = max{^ - 0}, (18) 

where j stands for the water level and is calculated according 
to the power constraint. 

As opposed to traditional wireless communication systems, 
AF relaying is not well suited for MI based transmissions, 
due to a high frequency-selectivity of the channel, especially 
if the relay is placed close to the source. This yields a narrow 
frequency band and a poor system performance. This behavior 
can be explained via Fig. [2] where normalized results for 
SNRm.rc(/) are shown using the optimal system parameters 
for the corresponding relay positions. If we neglect the passive 
relaying link due to much higher path loss, we can distinguish 
between three cases: 





















1) Relay is placed close to the source: Po n (f) is the 
dominant noise source, such that 

A 2 P S t(f)\H SD ,M )\ 2 


SNR 2 (/) = 


PDnU) 


(19) 


results. Due to an extremely frequency-selective path 
loss \Hsd, a(/)| 2 , the useful bandwidth is a few Hz. 


2) Relay is placed close to 

A 2 \H RR (f)\ 2 \H RD (f)\ 2 P Rn (f) 

noise source, such that 


the destination: 
is the dominant 


snr 2 (/) = 
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PRn(f) 
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results. Here, the frequency-selectivity is identical with 
that of the direct MI transmission, which provides the 
largest possible bandwidth of up to several 10 kHz. 

3) Relay is close to the middle: no noise component can 
be viewed as clearly dominant. The resulting bandwidth 
is larger than for case 1, but smaller than for case 2. 
Hence, the achievable data rate for AF relaying is maximum 
at the destination and the performance of this scheme is in 
general very poor. 


B. Filter-and-Forward relaying (FF) 

The FF relaying is similar to the AF relaying, however, 
a frequency-selective amplification coefficient A(f) is used 
instead of a constant A. As known from the literature (e.g. 
El. El), A(f) can be optimized in order to maximize 
the achievable data rate in wireless communication systems. 
Hence, the optimization problem can be formulated as follows: 

max [ log 2 (1 + SNR MKC (/)) d/, (21) 

A(f),p S t(f)m+ Jb 

s.t.: 1) P s = f Pst(f)df, 

Jb 

2) Pr= f Pm(f)df. 

Jb 

Unfortunately, it can be shown that this problem is not convex 
and its solution cannot be calculated using the well-known 
methods of convex optimization |fT6l . Hence, an iterative algo¬ 
rithm is proposed in Ell . In each iteration, two subproblems of 
© need to be solved. By solving the first subproblem, Pst(f) 
for a given relay transmit power spectrum density Pm(f) is 
optimized. Then, by solving the second subproblem, Pm{f) 
for a given power distribution Pst(f) is optimized. P R t(f) is 
related to the amplification A(f) by 

Pm(f) = \A(f)\ 2 (. P St (f ) \H SR (f)\ 2 + P Rn (f)) \H RR (f)\ 2 . 

( 22 ) 

According to El, both subproblems are convex, such that a 
closed-form solution for each subproblem can be found using 
a Lagrangian approach which is also given in fl4l . 


C. Decode-and-forward relaying (DF) 

In DF relaying, there is no linear mapping between the 
received signal and the transmitted signal at the relay, such 
that the two links (source —► relay) and (relay —» destination) 


can be decoupled and optimized separately. Due to a much 
higher attenuation of the passive relaying link, no diversity 
combining is needed at the destination node. The achievable 
data rate of the entire system can be given by 

Cdf,hd = -vain{C SR ,C RD }, (23) 

<24> 

w.- 4 + ^ ( g /)l > <25) 

The data rate ( |23l > is then maximized using the waterfilling 
rule applied to Pst.(f) and P R t(f ), respectively. 


D. Full duplex relaying mode (FD) 

In contrast to the traditional communication systems, the 
use of full duplex relaying mode is possible with MI based 
transmissions due to the mostly time-invariant transmission 
channels, especially in MI based WUSNs as mentioned in El- 
For this, the known signal transmitted by the relay is subtracted 
from the sum of transmitted and received signals at the load 
impedance R r of the relay. The received own signal at the 
load impedance of the relay during transmissions from the 
relay can be given by 

URr = Ur z* + (IM) 2 (iwL + AP rd ) ’ (26) 


such that the corresponding path loss is equivalent to 
\H RR (f)\ 2 from (HJ. Since U R and all system parameters are 
known to the transmitting relay, in principle it can perfectly 
separate the own transmitted stream from the other received 
signals. Of course, in case of sudden changes of the prop¬ 
agation characteristics (due to rainfall, etc.) M RR and M RR > 
may change, such that extraction of the desired signal becomes 
inefficient. However, this problem can be handled by a proper 
channel estimation, which is very efficient in WUSNs due to 
the stationary deployment. Hence, in the following we assume 
that a perfect signal extraction in full duplex mode is possible. 
Since the received signals in the same time slot from the relay 
and the source (via passive relaying) cannot be combined in 
case of FD mode, these signals impose additional interference 
to each other. We focus on the stronger signal, which is obvi¬ 
ously the signal coming from the relay to the destination. The 
second signal coming passively from the source is considered 
as interference. Hence, some equations shown in this section 
need to be changed accordingly: 

1) AF and FF relaying: For the AF and FF relaying 
schemes, we obtain 


SINR FD (f) 


1 + SNIl] 


SNR 2 (/) _ 

i , \A(f)\ 2 \Hsp.Af)\ 2 Pn„(f) \ ’ 
+ \HsnU)\ 2 PoM) ) 

(27) 


where SINR^^j(/) denotes the total signal-to-interference- 


plus-noise ratio in full duplex mode at frequency / as received 
by the destination. Here, the amplification coefficient A(f) is 










either frequency flat (AF relaying) or frequency-selective (FF 
relaying). With this, the achievable data rate can be given by 


C{AF,FF},FD — [ l°g 2 (l + SINR FD (/))d/. (28) 
J B 

2) DF relaying: For the DF relaying, only Crd needs 
to be adjusted, since Csr is not effected by the switch to 
the full duplex mode (assuming perfect extraction of the 
received signal is possible at the relay). Hence, we modify 
( |25) l similarly to ( |27| and ( |28j ): 


Crd,fd 




PRt(f)\H RD (f)\ 2 

PDn(f) + Pst(f)\H SD ,M )\ 2 



(29) 

Here, the signal transmission requires only one time slot. 
Hence, Cdf,fd = min {Csr,Crd,fd} holds. 


In general, the optimization algorithms shown in the previ¬ 
ous sections are valid for the full duplex mode as well. This is 
due to a very high path loss of the passive relaying link, such 
that the optimal system parameters and transmit filters do not 
change much compared to the half duplex mode. Therefore, 
we apply (|27j)-(|28j) only for the calculation of the final results 
after finding the optimal parameters Pst(f) and A for the 
AF and FF relaying. For DF relaying, ( f29| ) replaces ( |25[ i after 
the optimization of both Pst{f) and Pr 4 (/) as described in 
Section 1III-CI 

IV. Numerical Results 

In this section, we discuss the performance of different re¬ 
laying schemes based on the results obtained from a numerical 
evaluation of the proposed scheme. In our simulations, we 
assume a total transmit power budget of Ptotai = 10 mW. 
This includes the transmit power Pg at the source and Pr 
at the relay. The optimal values for Pg and Pr are found 
via full search under the constraint Pg + Pr = Ptotai- We 
utilize coils with wire radius 0.5 mm and coil radius a = 0.15 
m. The maximum number of coil windings N is 1000. The 
conductivity and permittivity of soil are a = 0.01 S/m and e = 
7eo for dry soil, respectively, where eo ~ 8.854 • 10~ 12 F/m. 
Since the permeability of soil is close to that of air, we use 
l_i = /i 0 with the magnetic constant fig = 4ir ■ 10“' H/m. For 
a reduced path loss, we assume that the axes of all coils are 
parallel to the direction of the signal transmission, thus, the 
polarization factors Jg r and Jrr from G are equal to two, 
which holds for horizontal axes deployment 0. In order to 
investigate reasonable positions of the relay, we choose the 
shortest distance between the relay and any other transceiver 
to be at least 3 m. 

At first, we show results for the achievable data rate using 
different relaying schemes in half duplex mode at different 
relay positions, see Fig. [3] Here, the performance for AF 
relaying differs from the conventional wireless relaying due to 
its maximum close to the destination. However, even placing 
the AF relay close to the destination does not provide a 
sufficient data rate to justify the use and deployment of this 
device, as we shall see later. The results for FF and DF relaying 



Fig. 3. Performance for different relay positions. 



correspond to our expectation, where the best performance is 
obtained if the relay is placed exactly in the middle. 
Furthermore, we show the optimal resonance frequency /o for 
different transmission distances and different relaying schemes 
in half duplex mode, see Fig. [4] For comparison, we also show 
the results for the direct MI transmission (without relay). As 
known from the previous works on magnetic induction based 
systems (e.g. 0 or El), the optimal resonance frequency 
decreases with increasing transmission distance. This behavior 
is confirmed here for the relay networks as well. Since the 
relay is placed in between the source and destination, the 
transmission distance of the longer link is smaller than the 
total transmission distance. Therefore, the optimal resonance 
frequency for all relaying schemes is higher than for the direct 
transmission. From the same reason, the resonance frequency 
for AF is lower than for FF and DF, which is due to the relay 
position in AF close to the destination. Moreover, the optimal 
resonance frequency for AF relaying converges to that of the 
























Fig. 5. Achievable data rate at the optimal relay position. 


direct MI transmission for increasing distance. For the relayed 
signal transmission, the curves obtained for the resonance 
frequency are not smooth over transmission distance, which 
is due to the finite precision of the optimization, where the 
resonance frequency is one of the many parameters jointly 
optimized via full search. 

The achievable data rate is shown versus transmission distance 
in Fig. [5] Obviously, AF relaying does not improve the per¬ 
formance of an MI based transmission, as mentioned earlier. 
Hence, this type of relaying should not be used in MI-WUSNs. 
On the other hand, we observe a very good performance 
of the FF relaying scheme, which almost reaches the upper 
bound given by the DF relaying. Hence, even non-regenerative 
relaying strategies can be useful for MI-WUSNs. However, DF 
is still 20% - 40% better than FF, which is due to the noise 
enhancement at the relay according to ( |22| > for FF. However, 
DF requires a much higher complexity of the signal processing 
and transceiver design, which is an important criterion for the 
design of MI-WUSNs. 

For full duplex relaying, the data rate is 50% - 65% larger than 
for the half duplex relaying in all considered relaying schemes, 
which makes this new technique very promising. Using DF in 
full duplex mode, data rate can be even increased by factor 10 - 
22 compared to the direct MI transmission with no active relay. 
In particular, even for distances up to 50 m, a data rate above 
1 Mbit/s can be achieved. Of course, in a practical system with 
finite modulation alphabet size and suboptimum filtering the 
resulting data rate may be somewhat lower. However, such an 
investigation is beyond the scope of this work. 

V. Conclusion 

In this paper, the potential of an active relaying technique in 
MI based WUSNs has been investigated, which comprises the 
traditional relaying of the data by means of amplification or 
decoding and the passive relaying inherent to MI based trans¬ 
missions. For this, a suitable system model has been proposed 
and the system parameters have been optimized in order to 


maximize the achievable rate. One of our main observations 
is, that the passively relayed signals have a very little impact 
on the result due to a much higher path loss. Secondly, a 
simple AF relaying performs even worse than the direct MI 
transmission without relay. This is due to a high frequency- 
selectivity of the path loss, such that the optimal position of 
the relay in this case is close to the destination. However, then 
the bottleneck link becomes dominant and the performance 
degrades. By using a frequency-selective amplification (FF 
relaying), the signal power can be redistributed in the relay, 
such that the frequency-selectivity of the channel can be 
partially compensated. Therefore, this relaying scheme shows 
a very promising performance, which is close to the upper 
bound given by the DF relaying. In addition, the possibility 
of utilizing a full duplex relaying mode with simultaneous 
signal transmission and reception has been discussed and the 
performance of this technique has been evaluated. 

References 

[1] I.F. Akyildiz, W. Su, Y. Sankarasubramaniam, and E. Cayirci, “Wireless 
sensor networks: A survey,” Comput. Netw. J., vol. 38, pp. 393-422, 
March 2002. 

[2] I.F. Akyildiz and E.P. Stuntebeck, “Wireless underground sensor net¬ 
works: Research challenges,” Ad Hoc Netw. J., vol. 4, pp. 669-686, 
July 2006. 

[3] I.F. Akyildiz, Z. Sun, and M.C. Vuran, “Signal propagation techniques 
for wireless underground communication networks,” Physical Commu¬ 
nication Journal (Elsevier), vol. 2, pp. 167- 183, September 2009. 

[4] Z. Sun and I.F. Akyildiz, “Magnetic induction communications for 
wireless underground sensor networks,” IEEE Trans, on Antennas and 
Propag., vol. 58, pp. 2426-2435, July 2010. 

[5] H. Jiang and Y. Wang, “Capacity performance of an inductively coupled 
near field communication system,” in Proc. of IEEE International 
Symposium of Antenna and Propagation Society, July 2008. 

[6] J.I. Agbinya and M. Mashipour, “Power equations and capacity per¬ 
formance of magnetic induction communication systems,” Wireless 
Personal Communications Journal, vol. 64, no. 4, pp. 831-845, 2012. 

[7] S. Kisseleff, W. Gerstacker, R. Schober, Z. Sun, and I.F. Akyildiz, 
“Channel capacity of magnetic induction based wireless underground 
sensor networks under practical constraints,” in Proc. of IEEE WCNC 
2013, April 2013. 

[8] S. Kisseleff, I.F. Akyildiz, and W. Gerstacker, “On Modulation for 
Magnetic Induction based Transmission in Wireless Underground Sensor 
Networks,” in Proc. of ICC 2014, May 2014. 

[9] E. Shamonina, V. A. Kalinin, K. H. Ringhofer, and L. Solymar, 
“Magneto-inductive waveguide,” Electronic Letters, vol. 38, no. 8, pp. 
371-373, 2002. 

[10] M. Masihpour and J.I. Agbinya, “Cooperative relay in Near Field 
Magnetic Induction: A new technology for embedded medical commu¬ 
nication systems,” in Proc. of IB2Com, December 2010, pp. 1-6. 

[11] S. Kisseleff, I.F. Akyildiz, and W. Gerstacker, “Interference Polarization 
in Magnetic Induction based Wireless Underground Sensor Networks,” 
in Proc. of IEEE PIMRC 2013 (SENSA Workshop), September 2013. 

[12] F.H.P. Fitzek and M.D. Katz, Cooperation in Wireless Networks: Prin¬ 
ciples and Applications. Springer Netherlands, 2006. 

[13] Y.W.P. Hong, W.J. Huang, and C.C.J. Kuo, Cooperative Communications 
and Networking: Technologies and System Design. Springer, 2010. 

[14] I. Hammerstrom and A. Wittneben, “On the Optimal Power Allocation 
for Nonregenerative OFDM Relay Links,” in Proc. of ICC 2006, June 
2006, pp. 4463^1468. 

[15] Z. Fang, Y. Hua, and J.C. Koshy, “Joint Source and Relay Optimization 
for a Non-Regenerative MIMO Relay,” in Proc. of IEEE Workshop on 
Sensor Array and Multichannel Processing, July 2006, pp. 239-243. 

[16] S. Boyd and L. Vandenberghe, Convex Optimization. Cambridge 
University Press, 2004. 

[17] S. Kisseleff, I.F. Akyildiz, and W. Gerstacker, “Transmitter-Side Channel 
Estimation in Magnetic Induction based Communication Systems,” in 
Proc. of IEEE BlackSeaCom 2014, May 2014. 








